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FIRST ORDER PERTURBATIONS OF AN ORBIT BY 
A MASS ANOMALY 

David E. Smith 
ABSTRACT 

The first order short period perturbations of a 
satellites position and the first order long period per- 
turbations of the orbital elements by a mass anomaly 
(mascon) are developed. 
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FIRST ORDER PERTURBATIONS OF AN ORBIT BY 
A MASS ANOMALY 


INTRODUCTION 

Mass concentrations on the moon have been discovered from the analysis of 
the orbits of lunar satellites (see, for example, references 1, 2, 3) and these 
mascons, as they are usually called, have made the determination of orbits of 
lunar satellites difficult and their prediction rather unreliable. 

The effect of a mass concentration is, in certain respects, similar to that 
of a third body except that the mass concentration is entirely inside the orbit of 
the spacecraft and is in a fixed relationship with respect to the primary body. 

This similarity, however, enables the perturbing effects of a mascon to be calcu- 
lated from a very restricted three body approach as might be employed, for 
example, in computing the first order effect of the moon on a near Earth satellite. 

This is the approach that has been followed here. The disturbing function of 
a mascon has been obtained and used to develop the first order short and long 
period perturbations of the orbit of a spacecraft. 

DISTURBING POTENTIAL OF A MASCON 

Let us consider the potential at the point P (x, y, z) due to a mascon, m, at 
(X, Y, Z) and another mass (M - m) at (xj, y x ,Zj ). Let the origin of the co- 
ordinate system be at the center of gravity of the two masses (see Figure 1) and 
let the distances of P from m, the origin and (M - m) be p, r, t respectively 
The potential (V) at P can therefore be written. 
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v = G (M - m) + G_m (1 

l P 

Further, because the two masses m and (M - m) are in a fixed relationship, the 
disturbing potential is the difference between the potential at P and the potential 
of the central force term alone at P. Now the potential of the central force term 
is 


CM 

r 


and hence the disturbing potential (A V ) can be written as 


A V 


= V - 


GM 


r 


= G (M -m) + Gm - GJrt (2) 

Pi P r 

From the evaluation so far undertaken the size of a mascon is of the order 
of 10 5 of the moon’s mass so we can justifiably neglect terms of order (nv'M) 2 . 

In addition, because the origin of the coordinate system is at the center of gravity, 
we have 


r i x i y i m f e x I 

— = — = — + terms ot order 

R X Y M 


m 


M 


(3) 


2 



and 


li< m (4) 

r M 

where r 1 is the distance of (M - m) from the origin. Therefore, neglecting terms 
of order (nFM) 2 we have from Figure 1 

— =-(l--cos<p) (5) 

Pi r r 


and equation 2 can be written 


AV =Gra 



l R cos <t> 

r _2 


(6) 


In addition, from Figure 1, we have 


cos 4> 


xX + yY + zZ 
rR 


(?) 


and therefore the disturbing potential of a mascon can be written 


AV = Gm 


1 

P 


1 

r 


(xX + yY + zZ) 


+ terms of order 



(8) 


3 



From the disturbing potential we can obtain the disturbing acceleration, of 
which the x-component is 


4“ (AV) = - G m 
ox 




(x -X) 


(xX + yY + zZ) 


(9) 


with similar expressions for the y and z components. In order to develop the 
perturbations to the elements we require the components of the disturbing 
acceleration with respect to the orbit. Let the component acting along the radius 
vector be R , the component acting perpendicularly to the orbit plane S and the 
component in the orbit plane acting perpendicularly to R be Q , then 


R = *1 ~ ( AV) + m, (AV) + rij (AV) 

ox dy dz 

s - ^ 2 |- (AV) + m 2 | (AV) + n, 7- (AV) 

dx dy OZ 

0 = T (AV) + m 3 | (AV) + n 3 |- (AV) 
J ax J °y oz 

where 


( 10 ) 


( 11 ) 

(12) 


(^ r rij ) , (t 2 . m 2 .n 2 ) and (^ m 3 n 3 ) 

are the direction cosines of K ,Q and S . 

From Figure 2 we have, following the method of Cook (reference 4), 


0 x y z 

■v x = — , mj =— , nj - — 

r r r 


( 13 ) 


-l~ = cos ft sin u - sin ft cos u cos i 
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m 2 = - sin 0 sin u + cos Q. cos u cos i 
n 2 = cos u sin i 

'i 3 - sin n sin i 
m 3 = - cos H sin i 

n 3 = cos i 

where 

ft is the right ascension of the node 
u is the argument of latitude 
i is the orbital inclination 


( 14 ) 


(15) 


If the position of the mascon is right ascension a and decimation 8 , then we 
also have 


X = R cos 8 cos a 
Y = 8 cos S sin a 

Z = R sin 8 

The distance p of the satellite from the mascon can be obtained in terms 
of r, R and 4> using the cosine rule which, upon expansion, can be written as 
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p 3 = ~3 + 3 7 C °S0~f (f) (l-5 COS 2 <p) _ 


( 17 ) 


With the aid of equations 9 through 17 we can write the major terms in the 
components of the disturbing acceleration as follows: 


R = — f— ^ [l - - (A 2 + B 2 )- 3 A B sin 2 u - - (A 2 - B 2 ) cos 2 u 
3 \2r / 2 v ' 2 v ‘ J 

t* ' ' w 


(18) 


= ^ AB cos 2u - I (A 2 - B 2 ) sin 2 U J 


(19) 


S _ Gm ^3R_ [A cos u + B sin u] 


(20) 


where 


A = cos § cos (a - fl) (21) 

B = sin i sin 8 + cos i cos 8 sin (a - fi) (22) 

C = cos i sin S - sin i cos § sin (a - fi) (23) 

The notation used above is the same as that used by Cook in reference 4 where 
details of the method employed here may be found. The terms neglected in 
equations 18, 19 and 20 correspond to terms of order (R/r) and above in equation 17. 
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PERTURBATION EQUATIONS 

The perturbations to the elliptic elements as a result of the disturbing ac- 
celeration are given by Lagrange's Planetary Equations. One form of these 
equations (references 4 and 5) accurate to order (m/M), is 


da 

dd 


2 r 2 
n 2 a p 


R e sin 6 + — 

r 



(24) 


de 

dd 



n 2 a 3 


R sin 6 + Q (cos 0 + cos E ) 


(25) 


dfi S r 3 sin u 
d6 n 2 a 3 p sin i 


(26) 


di_Sr 3 cosu (27) 

d9 n 2 a 3 p 


dec 

d 6 



- R cos 6 + (1 +— ) Q sin 6 
P 


dfi 

cos 


i 


de 


(28) 


da _ 2 r 3 R 

n 2 a 4 (1 - e 2 ) 1/2 


- (1 



da; 

_d6 


+ COS 



(29) 
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where 

a - M - j* ndt 

a is the semi-major axis 
e is the eccentricity 

u. is the argument of perigee 
M is the mean anomaly 
E is the eccentric anomaly 

r is the true anomaly 
n i s the mean mot ion 

and 

p = a (1 - e 2 ) 

In addition to the Keplerian elements it is of interest to know the perturba- 
tions to radial distance (r), position in the orbit and direction of travel. The 
latter is effectively given by the short-period perturbations of i for mascons 

near the equator but the perturbation equations for r and position in the orbit 
need to be derived. 

The radial distance is given by 


P 

r = 

1 + e cos 8 


( 30 ) 
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and hence 


dr 

d^ 


er 

sin 

P 


r 

a 


da 

dr 


r -v de 

- — (2ae + r cos r) — ; 

p do 


(31) 


Substituting for da/dr and de/d 6 from equations 24 and 25 and simplifying leads 
to 


dr er 2 sinr _ r 4 
dr p n 2 a 3 p 


sin 6 cos 0.R + ^2 +— 


cos 


ae 


p cos r 
re 



(32) 


for the instantaneous perturbations of radial distance. 
The mean anomaly M has been defined as 


M = a + 



(33) 


and hence 


dM do | dn . 

d 6 d 6 J d& 


do 

dd 


+ 3 I 2 da 

+ 2a 3 (1-e 2 ) 1/2 ^ 60 


dd 


dtr 

dd 


WzjVl f r 4 (R e sin 8 +- Q) d0 

n 2 a 3 p 2 J 


(34) 
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Substituting for da/d ~ from equation 29 leads to 


^ + <1 -eV /! 
da 


da, . d Cl 

+ COS 1 

_d& d& _ 


2 r 3 


(1 - e 2 ) 


2 \ 1/2 


n 2 a 3 p 


- 3(l-e 2 ) 1/2 

S --7 

n 2 a 3 p 2 


Jr 4 e sin 9 + Q j 


(35) 


where the left hand side is effectively the perturbation to the along track 
position. 


SHORT PERIOD PERTURBATIONS 

Substituting for R and Q in equation (32) for the perturbations to the radial 
distance and integrating leads to 


8 r 




cos 26 


-|a B s i n 2 u + -i (A 2 - B 2 ) cos 2u^- 


+ |AB sin 2 a; + — (A 2 - B 2 ) cos 2a;^^ + ^ cos 4 6 + cos 3 6 - cos 5 dj 

+ -|a B cos 2 a;-— (A 2 - B 2 ) sin 2a;j- ^ sin 40 + sin 3 6 - ~ sin 5 &J 


+ terms of order e 2 


t 


(36) 
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where 8 r is the radial perturbation of r. 

Similarly, substituting for R and Q in equation 35 and integrating twice, 
where necessary, leads to 


S M + (1 _ e 2 ) 1/2 


8 cc + cos i 8 Q 


= - 3 


(5) {ff < * 2>,/2 {> -! <A2 + b2 >} (s - r sm e) 


+~|a B cos 2u (A 2 - B 2 ) sin 2 u| 


e^AB sin 2 a. + -~(A 2 -B 2 ) cos 2a.'|-|^ sin 6 - sin 3 <9 


+ e |a B cos 2 a; - ~ (A 2 _ B 2 ) sin 2 a;^ cos 6 + ^ cos 3 


(37) 


The change in direction of travel of a spacecraft can be represented by the 
change in orbital inclination. Substituting for S in equation 27 and integrating 


8i = c 


2 \M / \P 


A 6 + — (A sin 2 u - B L os 2 u) 


2 

+ e sin 6 (1 - — sin 2 &) (A cos 2 a; + B sin 2 a) 


+ — e cos 3 6(A sin 2 a, - B cos 2 a;) 


(38) 
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Equations 36, 37 and 38 give the first order perturbations of the radial and 
along track positions and the direction of travel of a satellite. The maximum 
numerical size of the short periodic terms (those containing u or a) for a near 
lunar satellite can be seen to be of the following order for a mascon of 10" 5 of 
the moon's mass: 


b r 



~15 meters 


5M + occ + cos i Sfi 


~4 arc seconds 


Si 


3 

8 



arc seconds 


LONG PERIOD AND SECULAR PERTURBATIONS 

Substituting for R, S and Q in equations 24, 25, 26 and 28 and integrating from 
0 to 2 77 with respect to 6 , leads to 


Aa = 0 


(39) 


Ae = 0 


AD = 3W-M-YBC/sin i 


mV Rj 


Aa ; + cos i AG - - 2m ( — p ^ 


M/\py 


1 - ~ (A 2 + B 2 ) 


(40) 

(41) 

(42) 
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where the A quantities represent the change in the parameters per revolution of 
the satellite. 

1 he long period perturbation of the inclination can be obtained from equation 
38 by neglecting terms in sine and cosine <5 and replacing 6 by 2 n. We there- 
fore have 


Ai 



(43) 


Similarly, the long period and secular perturbations of the along track 
position can be obtained from equation 37. Hence 


AM + (1 


e 2 ) 1/2 


Ace + cos i Ai 1 


- 6r 



(1 - e 2 ) 1 ' 

P 


2 fl _ — 
2 


(A 2 + B 2 ) 


(44) 


Because there are no long period or secular perturbations of the semi- 
major axis (equation 39) we have 


A a = AM 

and hence equation 44 could have been obtained from equation 29 directly. 

For a close lunar satellite and (m/M) = 10' 5 the maximum magnitude of 
the first order long period and secular perturbations are 
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Afi-x,8 / sin i arc seconds'revolution 


Ao. + cos i Afl -v, 16 arc seconds / revolut ion 


Ai 8 arc seconds'revolut ion 


AM + Ao. + cos i A ^32 arc seconds/revolut ion 


CONCLUSIONS 

The first order perturbations of an orbit by a mass concentration have been 
developed and it has been shown that the maximum values of the secular and long 
period perturbations are about 300 and 80 meters in the along track and across 
track directions (1 arc second on the moon’s surface is about 10 meters) per 
revolution for a mascon of 10 5 of the moon's mass. In addition, short period 
perturbations have been found with maximum amplitudes of about 15 meters in 
height, 40 meters in along track position and 7 meters in the across track 
direction. 

The major assumptions that were made in the analysis were that the position 
of the mascon and the orbit of the spacecraft remained fixed in space over the 
integration period (one revolution). The former approximation can be justified 
because the moon rotates very slowly (about 0.5 degrees) per hour) and the 
second approximation is valid if the perturbations by all other forces are small, 
which is usually the case. 
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The main restriction on the whole theory is that it is only of first order. 
Neglecting Higher powers of (R/r) in the expansion of the disturbing acceleration 
limits the region of applicability of the theory. For example, if the mascon is 
near the surface of the body (R comparable to the radius of the body) then the 
errors in the acceleration are of order 10% for 10R. However, if the mascon 
is near the center of the main body, the theory is valid for even very close 
orbits and should indicate the true changes in the orbit. Qie further restriction 
is that the theory cannot be applied to circular orbits. The integrations performed 
were with respect to the true anomaly thus implying that a perigee could be 
defined and therefore the orbit was non-circular. 

In addition to the application of the equations that have been derived to 
certain orbital computation problems, the analysis enables the approximate 
magnitude and character of mascon perturbations to be easily assessed, which 
could be useful in preliminary studies of orbit problems. 
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Figure 1. Coordinate System for Point Mass 
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Figure 2. Satellite Orbit, Mascon and Components of Disturbing Acceleration 
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